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ABSTRACT 


The Spectral Analysis Program (SAP) was developed to provide the Manned 
Spacecraft Center with the capability of computing the power spectrum of a 
phase or frequency modulated high frequency carrier wave. Previous power 
spectrum computational techniques have been restricted to relatively simple 
modulating signals because of excessive computational time, even on a high 
speed digital computer. The present technique uses the recently developed 
extended fast Fourier transform and represents a generalized approach for 
simple and complex modulating signals. The present technique is especially 
convenient for implementation of a variety of low-pass filters for the 
modulating signal and bandpass filters for the modulated signal. 

The Spectral Analysis Program was written in Fortran IV for the 
UNIVAC 1230 digital computer. 
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1. INTRODUCTION 


The purpose of this guide is to provide the user with the necessary 
information to 1) understand the general computational approach used in the 
Spectral Analysis Program (SAP), 2) set-up the input parameters in the appro- 
priate tape or card format, and 3) execute the program. The detailed mathe- 
matical formulation is presented in Reference 1 (Final Project Report). 

A programmer's manual, containing the software flow diagrams and listings, 
is presented in Reference 2 (Spectral Analysis Program, Volume II - Pro- 
grammer's Manual). 

This guide has been written in two basic sections - the Program Descrip- 
tion section and the User's Information section. The Program Description 
section summarizes the mathematical- model and presents the basic software 
structure in terms of a simplified block diagram. A discussion of limiting 
cases, accuracy, and operational constraints is also given in this section. 

The User's Information section describes 1) the magnetic tape and card input 
required to successfully set-up and execute the program and 2) the output 
power spectrum print and/or plot. A sample case and timing data are included. 

This guide also includes a description of the Plot Generation Program 
(PLTGEN) , used in conjunction with SAP, which generates an input tape for the 
EAI Dataplotter. 



2. PROGRAM DESCRIPTION 


2.1 PROGRAM DEFINITION 

The Spectral Analysis Program (SAP) computes the power spectrum of an 
angle modulated high frequency carrier wave. The program is designed to 
1) provide capability for performing spectral analyses of a wide variety of 
angle modulated signals likely to be encountered in the testing and evaluation 
of angle modulated communication systems, 2) be operationally compatible’ wi th 
the Electronic Systems Test Laboratory (ESTL) UNIVAC 1230 digital computer, 
and 3) generate a power spectrum data tape compatible with the Plot Genera- 
tion Program (PLTGEN). PLTGEN generates the appropriate commands to drive 
the EAI 3440 Dataplotter and create the power spectrum plot. PLTGEN and the 
associated operating instructions are presented in Section 4. 

A simplified block diagram of SAP is presented in Figure 2.1. SAP 
either accepts an externally supplied modulating signal existing on magnetic 
tape or internally generates a user specified modulating test signal. The 
user specified signal may be selected from 1) a sum of sinusoids containing 
up to 25 frequency components (the amplitude, phase, and frequency of each 
component are user specified), 2) a periodic' four level signal in which the 
incremental amplitude levels and switching times are specified, and 3) a 
periodic square wave in which the amplitude and period are specified. After 
acceptance or generation, the input modulating signal may be low-pass filtered 
using either a Butterworth, Chebyshev, or user defined filter. For the 
Butterworth and Chebyshev filters, the user selects the filter order, cut- 
off frequency, and the amount of ripple (Chebyshev filter only). For the 
user defined filter, the filter transfer function is defined by the linear 
interpolation of the magnitude and phase for up to 50 frequency break-points 
supplied by the user. After low-pass filtering, the modulating signal is 
multiplied by the modulation parameter and then exponentiated to simulate 
the exponential modulation process. Next, the modulated signal may be band- 
pass filtered, using the same filters mentioned above. In addition, the 
bandpass center frequency may be specified. Finally, the power spectrum is 
computed, two output tapes are generated, and the output spectrum is printed. 
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As stated above, two output tapes are available at the conclusion of 
a SAP execution. The first tape is called the SAP plot tape and is used as 
input to PLTGEN. PLTGEN, in turn, generates another tape which is mounted 
on the EAI plotter unit and drives the plotting pen. The second tape is 
called the SAP spectrum tape and contains the Fourier transform of the 
modulated signal. This tape is re-input to SAP after the first SAP execu- . 
tion when displays of additional spectral regions are desired, as discussed 
below. 

The computed power spectrum may contain up to 16,384 spectral lines at 

a frequency resolution Af= yg ^§ ' 4 A ' t > where At is the sampling interval in 

seconds. If desired, a specified number of adjacent spectral lines of the 
computed spectrum may be averaged to provide a power spectrum of lower 
resolution. 

SAP is implemented based on the modular software design concept. The 
modules closely resemble the physical process being simulated. The input 
signal for a given module resides on one of two magnetic tapes; the output 
signal is computed and stored on the second tape. The second tape then 
contains the input signal for the next module. As a result of this signal 
processing technique, SAP may be restarted at an entry point corresponding 
to the last, completely executed module. 

Several power spectrum print and plot display options are available to 
provide maximum display utility to the user. Each point or plot may contain 
up to 480 spectral lines (corresponding to a maximum plot resolution of 20 
lines per inch). In the carrier centered display, the carrier is centered 
along the frequency axis and the spectral region associated with 240 lines 
to the left of the carrier and 240 lines to the right of the carrier is 
displayed. In other displays, the starting frequency is user specified and 
the spectral region associated with the starting frequency plus the next 
479 lines is displayed. Spectral regions larger than 480 lines may be 
displayed by additional plots in which the starting frequency is the last 
frequency of the previous plot. The entire SAP program need not be re- 
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executed in this case; only the last part of the program (from entry point 
number 10) is re-executed. Thirty-five plots are required to display the 
maximum size computed spectrum of 16,384 lines. 

All parameters describing the plot are input to SAP; PLTGEN does not 
require any user controlled input parameters. 

2.2 METHOD 

2.2.1 Mathematical Formulation 

The modulated signal e{t) generated by the process of angle modulation 
may be written 

e(t) = cos[w c t + 8<|>(t)] (1) 

where 4 >(t) represents a time varying phase function g(t) normalized to unit 
peak amplitude, i .e. , 


. fx.\ 

nw 


g(t) 



< 1 . 


The constant B is termed the modulation parameter and corresponds to the 
maximum phase deviation of the carrier, u is the carrier angular frequency 
and a) t + B<f>(t) is the angle of the modulated signal. 

For phase modulation (PM), <f>(t) is the normalized modulating signal and 
the peak phase deviation is represented by the modulation index 8. For, 
frequency modulation (FM), ( t ) is the integral of the normalized modulating 

signal. The instantaneous frequency deviation from the carrier frequency is 

f i ' f c = It * ,(t) = 

where <}>'(t) is the time derivative of <j>(t) and |<f>'(t)| <_ 1. Thus, for FM, 

B is the peak angular frequency deviation of the carrier. 


2-4 


By writing e(t) in exponential form as 



the components of e(t) identified with positive and negative frequencies 
can be associated with the positive and negative exponential terms. The 
spectrum of e(t) is obtained by taking its Fourier transform. 

E(u>) = E + (oj) + E_(u)) = F[e(t)] (3) 

= 1 iFEe^c 1 ♦ Ffe-^c* ] . 

But, 

F [^“c* e jB4>(t) ] = 6 (a) - a) c ) * F[e jei|>(t) ] , (4) 

where * indicates convolution. Consequently, 

E(cu) - \ «U - » c ) * F[e js * (t) ] 

*\ S(« + « c ) * F[e‘ je * (t >] .. (5) 

For a typical situation encountered in a communications system, the carrier 
frequency, ui Qt is large with respect to the spectral extent of F[e J ' B ^ t ^]. 

In such a situation, the spectra associated with the two exponential terms 
do not overlap. The sketch of Figure 2-2 illustrates a typical situation. 

Since the two spectra are identical except for frequency location, 

computation of the spectrum for only the positive exponential term in 

Equation (2) is sufficient. In addition, the multiplicative term, e Ja) c t J 

serves only to position the spectrum of F[e JB<t> ^] in region around u . 

c 

Consequently, the basic analytical problem to be resolved is the computation 
of the quantity 

V + M = F[e jM(t) ]. (6) 
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Figure 2-2. Typical Spectrum of an Angle Modulated Signal 


The desired output spectrum, E(w), is then 
E(oj) = j V + ( a) - w c ) + -i V + ( u) + u c ) 


(7) 


For those instances where the separation of the positive negative frequency 
terms of E(w) is not distinct, the two separate spectral terms can be combined 
by computing 


E(o)) - y V + (w - oj c ) + V_(w + u c ) , 


(8) 


where V_(u) = F[e - ^ ( t ) j ^ 
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.The previous discussion presented the overall mathematical description 
of SAP. The specific mathematical description for each module or subroutine 
is discussed below. 

a) TSGEN - Computes the internally generated modulating signal 
4>(t). These signals are: 

i ) Sum of sinusoids 

M 

4 >(kAt) = £ A. sin (2irf .kAt + a) 
j=l J J J 

where 

x i_ 

4 > ( kA t ) = k tn sample of signal 
At = sampling interval 
M = number of sinusoids 
A, = amplitude of j*"* 1 sinusoid 

J 

XL 

f . = frequency of j sinusoid 

v 

XL 

0 . = phase of j sinusoid 

J 

ii) Periodic four level modulating signal (illustrated 
on Figure 2-3) 

z 

<j>(kAt) = £ A. 

j=i j 

where 

A. = change in signal at t. time point 

J J 

XL 

t. = time point of j 1, change in the signal 

J 

£ = index defined by k. At, and the t.'s 

J 

such that t £ _.j kAt - mt g _< T where 

tjj = period of signal, m = 1, 2 

T = record length, and t Q e 0. 
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SIGNAL REPEATS 
EVERY t g SECONDS 


Figure 2-3. Periodic Four Level Test Signal 
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Figure 2-4. Periodic Square Wave 



iii) Periodic square wave modulating signal (illustrated 
on Figure 2-4) 

4>(kAt) = + A-j 

where the + and - signs are interchanged every 
t-j seconds, i.e., period of square wave is 2*t^ . 

b) ISAR - resamples (using linear interpolation) user supplied 
modulating signal and generates a signal compatible with the 
program 

<W (kAt > = *in (KAT > 

+ U in (KAT + aT) - 0. n (KAT))*(kAt - KAT)/AT 

where <J>.j n (KAT) = the K*'* 1 sample of the input signal 

AT = sampling interval for the input 
signal 

X L 

4> 0U t(kAt) = the k sample of the output signal 

At = sampling interval for the output 
signal 

(=AT/KSR) 

and the relation between k and K is 
KAT < kAt £ (K+1 )aT 

c) EFFT - computes the discrete Fourier Transform of inverse of 
a signal. The transform and inverse are defined by the pair: 

1 N- 1 ik 

x(j) = TT E x(k) V K 

N k=0 N 

N-l ik 

x(k) = E X(j) wj k 
j=0 N 
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for j = 0 , 1 , . . . , 
defined by: 


N - 1 and k = 0, 1 


N - 1 . i s 


W N = exp(2iri/N) 
where i = f-Y . 

The details of this calculation are given in Appendix A. 

d) FILTER - filters the modulating and/or modulated signal by 
multiplying the transform of the signal by the transfer function 
of the appropriate filter. The equation is 

W f > ■ H < f > S in< f > 

where 

S. (f) = Fourier Transform of input signal 

H (f ) = Fourier Transform of the impulse response 
of the filter (i.e., the transfer function) 

S out (f) = Fourier Transform of the output signal 

f = frequency 

e) TRFN - routine which evaluates the transfer function H(f) in 
the FILTER module. Three possible transfer functions are 
available - a user defined filter, a Butterworth filter, or 
a Chebychev filter. A detailed description of TRFN is given 
in Appendix B. 

f) MOD - computes the modulated signal e(t) 

e(t) = exp[-je<j>(t)] 
exp = exponential 
j = /=T 
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8 = modulation parameter 


<j>(t) = modulating signal for phase modulation 
or integral of modulating signal for 
frequency modulation 


g) PLOT - generates power spectrum data for display by: 

NAVE 

P'(kAf) = aE P(kAf + iAf - Af) 
i=l 

where 

Af = frequency increment 

a = scale factor computed by [NAVE*MAX(P(kAf ) )] 

P(kAf) = power in signal at frequency (kAf) 

P'(kAf) = average of P over NAVE frequency 

increments beginning at frequency kAf 

NAVE = number of frequency increments over 
which output data is averaged 


2.2.2 Software 


Cementation 


The Spectral Analysis Program (SAP) software is implemented using the 
modular design concept. Each module performs a specific operation on the 
ipput signal to yield an output signal. These operations correspond closely 
to the physical process being simulated. 


The overall program structure is presented by Figure 2.5. All modules 
and subroutines are shown. These routines and the associated functions 
were defined previously. 


SAP uses two tape units alternately for the input and output of the 
software modules. Each module accepts the signal stored on the input tape 
(termed the data tape), processes the signal, and stores the output signal 
on the output tape (termed the spare tape). The output tape then functions 
as the input tape for the next module, i.e., the program switches the tape 
unit labels for both tapes. 
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Key to Subroutine 
Names and Functions 


SAP 

TSGEN - 

ISAR - 

TRANS - 
EFFT - 

FILTER - 
TRFN - 


Driver for the 
Spectral Analysis 
Program 

Test Signal Gener 
ation 

Input Signal and 
Resampling 

Transpose 

Extended Fast 
Fourier Transform 

Filtering 

Transfer Function 
Calculation 


MOD - Modulation Operation 

PLOT - Plot and Output 

Generation 



Figure 2-5. Spectral Analysis Program Functional Block Diagram 





















The points in the driver program at which tape unit label switching 
occurs are called entry points. Ten entry points exist and the program 
locations are shown in Figure 2.5. The entry points are labeled 100, 200, 
300, 400, 500, 600, 700, 800, 900, and 1000 on the output listing. As each 
entry point is passed during program execution, the following messages are 
printed: 

PROGRAM HAS PASSED ENTRY POINT XXXX . 

DATA IS CURRENTLY ON UNIT XL 

SPARE TAPE IS CURRENTLY ON UNIT X2. 

THERE ARE X ROWS AND Y COLUMNS OF DATA. 

EXECUTION TIME = X MINUTES, Y SECONDS. 


The modular software design provides a convenient recovery capability 
should a machine or other failure occur. To recover after a failure, the 
user specifies the last entry point passed, the data unit number, and spare 
tape unit number as provided by the last set of entry point messages. The 
driver program then jumps to the specified entry point and execution of SAP 
continues. (Note: the entry points are labeled 100, 200, 300, 400, 500, 

600, 700, 800, 900, and 1000 on the output listing; to re-execute SAP, the 
input entry points are 1, 2, 3, 4, 5, 6, 7, 8, 9, and 10, respectively.) 

2.3 LIMITATIONS 

This section lists SAP operational constraints and discusses some 
special considerations relative to the analytical formulation of a problem. 

2.3.1 Operational Constraints 

SAP operational constraints are listed below. 

. The total number of data samples N = NCIDT*NRIDT must not exceed 
16,384. NCIDT is the number of columns of data and NRIDT is 
the number of rows of data. 
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. The total number of data samples N = NCIDT*NRIDT must be an 
integral power of two; i.e., N = 2 n , where n = 1, 2, 3 

13, 14. 

. The maximum value of NCIDT or NRIDT is 256. (For example, 
if N = 16,384, the following factors are acceptable: 

N = NCIDT*NRIDT = 256*64, 128*128, and 64*256. If N = 4096, 

N = 256*16, 128*32, 64*64, 32*128, and 16*256.) 

. The sum of sinusoids modulating signal generated internally 
' must contain 25 or fewer sine wave components; i.e., M £ 25. 

. Use of the Butterworth filter is restricted to filter orders 
one through six; i.e., N0RD1 or N0RD2 = 1, 2, 3, 4, 5 and 6. 

N0RD1 and N0RD2 are the filter orders of the low-pass and 
bandpass filters, respectively. 

. Use of the Chebyshev filter is restricted to filter orders 
one, three, and five; i.e., N0RD1 or N0RD2 =1,3, and 5. 

Even order Chebyshev filters are generally physically 
unrealizable for equal source and load impedances. 

. Use of the user filter is restricted to 50 break-points; i.e., 

NF1 or NF2 _< 50. NF1 and NF2 are the number of frequency 
break-points for the low-pass and bandpass filters, respectively. 

2.3.2 Special Considerations 

Experience indicates that the user should be familiar with the leakage 
effect, aliasing, and the Gibbs phenomenon. These considerations are briefly 
summarized below. Refer to the Final Project Report (Reference 1) for 
additional discussion and references. 

Leakage Effect . SAP is based on the use of the discrete approximation 
to the continuous Fourier transform. The discrete Fourier transform (DFT) 
approximation is exact only when the input signal is periodic and the fre- 
quency components are harmonically related. Otherwise, a leakage effect 
exists. This leakage effect exists because a finite record is formed con- 
taining only T seconds of the actual signal x(t). Forming a finite record 
may be viewed as multiplying the signal by a rectangular window function 
g(t) with unit amplitude during the T second interval and zero magnitude 
at other times. Thus, 
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and 


x 1 (t) = g(t) x ( t) 


x ' ( t) = x(t) if ^ <_ t <_ j 
=0 if | t 1 

Multiplication in the time domain corresponds to convolution in the transform 
domain; i.e., 


X 1 (w) = G(oi) * X(w). 

Since G(w) is the Fourier transform of g ( t ) and has the form (sin cj)/io, the 
convolution results in a blurring of frequency components. For example, 
if x(t) is an infinitely long cosine waveform of frequency ia Q> the Fourier 
transform yields two impulses that are symmetric about zero frequency 
(Figure 2-6a). If, however, a finite length record is transformed, these two 
impulses will be convolved with G(w) = (sin w)/ u, yielding the results 
presented in Figure 2-6c. Thus, some of the energy of the impulse "leaks" 
through the sidelobes into adjacent frequencies. 

Aliasing . The term "aliasing" refers to the fact that high-frequency 
components of a time function will appear as low frequencies if the sampling 
rate is too slow. If the sampling frequency is f = 1/At (where At is the 
sampling interval in seconds), a "folding frequency" exists defined by 
f^ = -f $ /2 . An input signal x(t) will be represented correctly if the 
highest frequency component is less than the folding frequency. A frequency 
component 5 Hz higher than the folding frequency will appear as a frequency 
5 Hz lower than the folding frequency. Aliasing may be avoided by sampling 
the signal at a rate at least twice as high as the highest frequency present; 
i.e., at the Nyquist rate. 
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Aliasing is of particular significance when the input signal is not 
band limited. For example, a square wave contains an infinite number of 
frequency components and can not be sampled at the Nyquist rate. The ampli- 
tudes of the components, however, decrease with increasing frequency. Thus, 
the number of components with significant (in terms of user requirements) 
magnitudes is limited to a finite number. Aliasing is always present when 
sampling a square wave, but may be reduced in significance by increasing the 
sampling rate. The sampling rate, however, may not be made arbitrarily large 
because the SAP capability is limited to 16,384 data samples. 

Gibbs Phenomenon" . The Gibbs phenomenon exists when a Fourier series 
is used to represent a function which has a finite number of discontinuities. 
This phenomenon has been experienced during SAP executions when the modulating 
signal was a square wave. In this case, the square wave shown in Figure 2-7a 
is approximated by the wave shown in Figure 2-7b and the power spectrum is 
computed accordingly. 

As shown by the latter figure, ripples exists in the vicinity of the 
discontinuities. According to Fourier theory, the trigonometric series 
approximation of the square wave improves in a mean square error sense as 
the number of terms n in the series increases. The theory, however, also 
indicates that the ripples exists even for very large values of n. As n is 
increased further, the ripples show a proportionately increased rate of 
oscillation, but the ripple amplitudes relative to the magnitude of the 
discontinuity remains the same. As n -* », the ripples are compressed into 
a single vertical line at the point of discontinuity. The square wave can 
never be approximated in the vicinity of the discontinuity with an error 
less than 18 percent. 

2.4 ACCURACY 

The Spectral Analysis Program was verified at the individual module 
level and at the overall software package level. Verification was accom- 
plished using test cases with known theoretical results... For example. 
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<f> (t) 


— t 


(a) True Square Wave Modulating Signal 


♦ (-t) 



(b) Approximate Square Wave Modulating Signal 


Figure 2-7. 


Gibbs Phenomenon 


2 - 



at the overall software package level, the computed power spectrum was com- 
pared with theoretical results for phase and frequency modulating signals 
consisting of sine waves and periodic square waves. Results indicated a 
minimum of six^ significant figures accuracy for the power spectrum expressed 
in decibels. 

This high degree of accuracy is a direct result of using the extended 
fast Fourier, transform (EFFT) technique to compute Fourier transforms. . 

Use of the EFFT for complex modulating signals is advantageous because 1) the 
number of multiply-add operations is reduced very significantly compared to 
non-FFT power spectrum computational techniques and 2) the EFFT accumulates 
a long series of floating point numbers optimally. The optimal accumulation 
method is to add pairs of numbers, then pairs of the resulting sums, etc. 

The error generated by the EFFT was investigated experimentally by performing 
the EFFT of a unit impulse, performing the inverse EFFT, and then examining 
the difference between the original and reconstructed sets of data. At 
least seven significant figures of agreement, depending on the array size, 
was obtained. 
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3. USER'S INFORMATION 


3.1 INPUT DESCRIPTION | 

I 

The subsections below define the input parameters for the Spectral i 

Analysis Program, the formats for the inputs, the overall card deck setup, 
and the necessary detailed machine operating procedures. 

3.1.1 Input Parameters i 

For convenience, the inputs for the SAP are divided into a number of 
categories. The categories are general inputs, inputs to define the test 
signal or the input signal, inputs to define filtering for both the modu- 
lating and modulated signals, and inputs to define the plot parameters. The 
input symbols and definitions are given by Table 3-1. 

3.1.2 Formats for Inputs 

Input to SAP consists of a set of data cards and a magnetic tape (if 
ISOPT = 0) containing an externally generated modulating signal. The card 
input parameters, card group description, and format are described in 
Table 3-2. (A card group is a set of one or more cards that relate to a 
given computational function. For example, if the modulated signal is to 
be filtered with the user filter, the bandpass user filter group of cards 
must be included. If not, this group of cards is omitted.) 

Each SAP execution must contain the first three card groups in Table 
3-2; i.e., each execution must contain a title card (which may be blank), 
two integer data cards, and a fixed point data card. If 1) an externally 
generated modulating signal is to be used and 2) neither the modulating 
signal nor the modulated signal is to be filtered, additional card input is 
not included; otherwise, one or more card groups are required. Card group 
requirements are presented in Table 3-3. 

. The number of cards required for a 'group is given by the number of 
formats given for that group in Table 3-2. For example, two formats are 
given for the integer data card group. Thus, two integer data cards are 
required. The first card contains 16 parameters, each occupying five card 
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Table 3-1. Inputs for Spectral Analysis Program 


General 

IEP 


IDT 

1ST 


NCIDT 

NRIDT 


ISOPT 


IFM 


IFOPT 


- defines entry point to the driver by 
'1 - nominal run 

z: m 

K - restart entry point (as shown in Figure 2-1) 

- defines data tape and spare tape units. For normal 
cases, IDT = 1 and 1ST = 2. For restart capability 
(IEP f 1), signal data resides on IDT. 

- defines the number of points, N, in the test signal 
by N = NCIDT*NRIDT. Each must be a power of 2. For 
efficiency, NCIDT and NRIDT should be approximately 
equal (for example NCIDT = 64 and NRIDT = 128). For 
restarting (IEP f 1), NCIDT and NRIDT define the num 
ber of columns and rows, respectively, of the arrays 
on tape IDT. 

- defines signal option by 

'O, user input signal on tape IDT 

1, sinusoidal test signal 

n - 

2, periodic four level gray test signal 
. [3, periodic square wave test signal 


- defines modulation type by 



1 , FM modulation 

2, phase modulation 


- defines filtering options 

0, no filtering 

1, filter modulating signal only 

=5 « 

2, filter.modulated signal only 

.3, filter modulating and modulated signal 
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Table 3-1. Inputs for Spectral Analysis Program (Continued) 


BETA 


DELT 

FC 

IZPRT(I) 


defines modulation index by 

if IFM = 1, BETA = maximum frequency deviation (Hertz) 

if IFM = 2, BETA = maximum phase deviation (radians) 

time separation (in seconds) between points in test 
signal or between user input signal prior to resampling 

carrier frequency (in Hertz) 

special print control 

0, no special print 

1, special print at indicated entry point 


I 


ENTRY POINT 


1 

2 

3 

4 

5 

6 


2 

4 

5 

6 
8 
9 


Test Signals 

a) Sinusoids (ISOPT = 1) 

M - number of sinusoids (_< 25) 

A(I) - amplitude of 1^ sinusoid 

F(I) - frequency (in Hertz) of I th sinusoid 

THETA(I) - phase angle (in degrees) of I th sinusoid 

For IFM = 1, A(I) is the peak frequency deviation associated 
with each sinusoid 

For IFM = 2, A(I) is the peak phase deviation associated 
with each sinusoid 

Note: If ISOPT = 1, the input specification for BETA 

must be omitted 
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Table 3-1. Inputs for Spectral Analysis Program (Continued) 


b) Periodic four level gray (ISOPT =2) 

TP ( I ) - I = 1, 2, 3, 4 - defines the time points (in seconds) 

of the breaks or changes in the signal 

A(I) - I = 1, 2, 3, 4 - defines the corresponding changes in 
the amplitude of the signal 

TP ( 5 ) - defines period of signal 

c) Periodic Square Wave (ISOPT = 3) 

TP ( 1 ) - period (in seconds) of the square wave 

A(l) - amplitude of square wave 


Input Signal 


NR1 - number of words per record on input data tape 

MPT nnmkAM A-f iAAAAV»rlr aa mM if A +• fan/5 

no i ~ I iuiiiw 1 1 oi i cowi uo vii inivuv uubu 


KSR - defines sampling rate for resampled signal by DELT/KSR 


Filtering (depending on IFOPT - zero, one or two filters need to be 
specified) 

a) for the modulating signal 

IFTYP1 - defines filter type by 

1 , Butterworth 

=-2, Chebychev 

3, Input filter 
» 

N0RD1 - defines order of Butterworth of Chebychev filter 

(< 5 ) 

FCUT1 - defines cutoff frequency of Butterworth or Chebychev 
filter (in Hertz) 

RIP1 - defines ripple for Chebychev filter (refer to 
Appendix B) 
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Table 3-1. Inputs for Spectral Analysis Program (Continued) 


NF1 - for IFTYP1 = 3, defines number of points specified 
for input filter (<_ 50) 

TM AG 1(1) - for IFTYP1 = 3, defines magnitude of input filter 
transfer function at I^h frequency 

TPH1(I) - for IFTYP1 = 3, defines corresponding phase (in radians) 
at ith frequency 

FI (I) - for IFTYP1 = 3, defines frequency (in Hertz) at which 

magnitude and phase is specified 

- intermediate values of the magnitude and phase are 
generated by using linear interpolation over the 
input .values 

b) for the modulated signal 

IFTYP2 

N0RD2 

FCEN2 identical in definition as corresponding variables 

FCUT2 for modulating or test signal filter - but apply to 

RIP2 - filtering the modulated signal. The additional 

NF2 parameter FCEN2 defines the center frequency for the 

TMAG2(I ) bandpass filter. 

TPH2 ( I ) 

F2(I ) J 


Plot 

FSTRT - frequency (in Hertz) at which plot display begins; 

if < 0, then display with centered FC is generated. 

AMAX - scaling factor in output plot, if set to zero the plot 

module computes AMAX as maximum power component in 
total signal and scales data to be plotted by this value, 
if set positive then the input value is used to scale 
data for plotting. 

NAVE - number of points to be averaged for each output display 

point, each output point represents a bandwidth of 
NAVE*Af where Af = 1/(NRIDT*NCIDT*DELT) 
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Table 3-1. Inputs for Spectral Analysis Program (Continued) 


IPMODE 


IPLPOS 


ITITLE 


- controls printing and plotting by 
0, print and plot 

=' 1 , plot only 
2, print only 

- controls plot position on page by 
Oj top half 

~ « 

.1, bottom half 

- 30 characters for title on run and title on plot 
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Table 3-2. Card Input 


Card Group 


Description 

Parameters on Card 

Format 

Title 

(Up to 30 alpha-numeric characters) 

(30A1 ) 

Integer 

Data 

IDT, IEP, IFM, IFOPT, IFTYP1 , IFTYP2 , IPLPOS, 
IPMODE, ISOPT, 1ST, KSR, NCIDT, NCI, NF1 , 

NF2, NRIDT 

(1615) 


NR1 , NAVE, IZPRT(l) IZPRT(6) 

(815) 

Fixed Point 
Data 

AMAX, BETA, DELT, FC, FSTRT 

(5F10.0) 

Low Pass 
User Filter 

FI (1 ) , FI (2) , FI (3) FI (8) , 

FI (9) , FI (10) , Fl(ll), .... FI (16) , 
FI (17) , ..., Fl(NFl) 

(8F10.0) 

II 

II 


TMAGl(l), TMAG1 (2) , TMAG1 (3) TMAG1 (8) 

TMAG1 (9) , TMAGl(lO), TMAGl(ll), ..., TMAG1(16), 
TMAG1 (17) , ... , TMAGl(NFl) 

(8F10.0) 

II 

II 


TPHl(l), TPH1 (2) , TPH1 (3) , ..., TPH1(8). 
TPH1 (9) , TPHl(lO), TPH1 (11 ),..., TPH1(16), 
TPH1 (17) , ... , TPHl(NFl) 

(8F10.0) 

II 

II 

Low Pass 
Built-in 
Filter 

N0RD1 , FCUT1 , RIP! 

(15 ,2F1 0.0) 

Bandpass 
User Filter 

F2(l), F2(2) , F2(3) , .... F2(8), 

F2(9) , F2 ( 1 0) , F2(ll) F2(16), 

F2(17), ..., F2(NF2) 

(8F10.0) 

II 

II 


TMAG2(1), TMAG2(2) , TMAG2(3) TMAG2(8) , 

TMAG2(9) , TMAG2(10), TMAG2(11), ..., TMAG2(16), 
TMAG2(17) , ... , TMAG2(NF2) 

(8F10.0) 

II 

II 


TPH2(1 ) , TPH2(2) , TPH2(3), ..., TPH2(8), 
TPH2(9) , TPH2 (10) , TPH2(11), ..., TPH2(16), 
TPH2(1 7) , ..., TPH2(NF2) 

(8F10.0) 

II 

It 

Bandpass 

N0RD2, FCEN2, FCUT2, RIP2 

(I5.3F10.0) 


Built-in 

Filter 
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Table 3-2. Card Input (Continued) 


Card Group 
Description 

Parameter on Card 

Format 

Test Signal- 
Sinusoids 

M 

(12) 


A(l), A(2), A(3), .... A(8), 
A(9), A( 1 0) 9 A(ll ) , ...» A(16), 
A(17), , A (M ) 

(8F10.0) 

II 

II 


F ( 1 ) , F(2), F(3) , ...» F (8) , 

F ( 9 ) , F ( 1 0 ) , F(ll) 5 .... F(16) , 
F(17), ..., F(M) 

(8F10.0) 

II 

II 


THETA(l), THETA(2), THETA(3) , .... THETA(8) > 
THETA(9) , THETA(IO) , THETA (11 ) , .... THETA(16) , 
THETA(17) , .... THETA(M) 

(8F10.0) 

II 

II 

Test Signal- 
Multilevel 

TP ( 1 ) » TP ( 2 ) , TP(3), TP (4) , TP(5) 
A(1 ) , A ( 2) , A(3), A(4) 

(8F10.0) 

(8F10.0) 

Test Signal- 
Square 

TP, A 

(2F10.0) 
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Table 3-3. Card Group Requirements 


Card Group Requirement 


Title 

Always required 


Integer Data 

Always required 


Fixed Point Data 

Always Required 


Low Pass User Filter 

IFOPT = 1 or 

3, 

IFTYPl = 3 

Low Pass Built-in Filter 

IFOPT = 1 or 

3, 

IFTYP1 = 1 or 2 

Bandpass User Filter 

IFOPT = 2 or 

3, 

IFTYP2 = 3 

Bandpass Built-in Filter 

IFOPT = 2 or 

3, 

IFTYP2 = 1 or 2 

Test Signal-Sinusoids 

ISOPT = 1 



Test Signal -Multi level 

ISOPT = 2 



Test Signal -Square 

ISOPT = 3 
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columns (i.e., each parameter card field corresponds to five card columns). 
The second card contains eight parameters, each occupying five card columns. 
In this case, the first eight card fields must be used. The parameter must 
be key-punched on the card in the order shown by the table. 

The low pass user filter, bandpass user filter, and the test -signal- 
sinusoids card groups use a variable number of cards, as indicated by the 
" marks. For example, if the low pass user filter is to be used and NF1 = 8, 
then one card containing eight FI values, a second card containing eight 
TMAG1 values, and a third card containing eight TPH1 values are required. 

If NF1 = 50 (corresponding to the maximum allowable break-points), then the 
following cards inputs are required: six cards containing eight FI values 

plus one card containing two FI values, six cards containing eight TMAG1 
values plus one card containing two TMAG1 values, and six cards containing 
eight TPH1 values plus one card containing two TPH1 values. 

A magnetic tape containing an externally generated modulating signal 
is required if ISOPT = 0. This tape contains NCI number of records and NR1 
number of words per record. Words are written at 556 bpi in unformatted 
binary. 

3.1.3 Overall Deck Setup 

The overall deck setup and control cards required to execute SAP on 
the UNIVAC 1230 computer are shown in Figure 3-1. 

The UNIVAC 1230 control cards are described below: 

. FORTN 

This card specifies to the Compiler Executive that a FORTRAN 
program follows and that the program is to be compiled. FORTN 
starts in column one. 

. CONTROL ,T, BASE 646 

This card specifies that the desired mode of output to be used by 
the Compiling System is tape. A base of 646 is used for the main 
program. The base for subroutines is optional. -CONTROL starts in 
column seven. 
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Figure 3-1. Overall Deck Setup 


. MAIN PROGRAM NAME 

This card designates the main routine. MAINPROGRAM begins in 
column seven. 

. END 

This card specifies the end of the compilation. END begins in 
column one. 

. LOADT NAME 646 

This card specifies to the Compiler Executive that a previously 
compiled program, which has been output onto magnetic tape, is 
to be loaded into the computer starting at location 646. LOADT 
begins in column one. NAME begins in column eleven. "646" 
begins in column twenty-one. 

I 

. EXECUTE NAME ACCEPT 

This card indicates by name which program is to have control. The 
word EXECUTE begins in column one. The name of the program (NAME) 
begins in column eleven. The word ACCEPT is unnecessary unless 
it is desirable to operate the program under non-fatal error 
circumstances. ACCEPT starts in column twenty-one. 

. DATA 

This card specifies to the Compiler Executive that input data may 
follow, and that the program is to be executed. This card is 
necessary to cause the execution of the program even if there are 
no data cards present. DATA begins in column one. 

. ENDATA 

This card signifies the end of the data stream. END begins in 
column one. 

. SYSEND 

This card signals the end of the FORTRAN job and causes control 
to be transfered back to the Compiler Executive. SYSEND begins 
in column one. 

3.2 OUTPUT DESCRIPTION 

3.2.1 Printed Output 

The print output is a listing which includes 1) case title, 2) all 
input parameters, 3) tape unit residence for input and output tapes, module 
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execution times, and entry point passage, 4) frequency and power (in decibels) 
for 480 points of the computed spectrum. An example of the printed output 
is presented in Section 3.4. 

3.2.2 Plotted Output 

The plotted output is a plot of 480 points of the computed power 
spectrum. The plot ordinate is relative power expressed in decibels in the 
range 0 to -70 dB ; the abscissa is frequency expressed in Hertz, Kilohertz, 
or Megahertz. Each power component is represented by a vertical line 
extending from -70 dB to the component value. The case title, carrier fre- 
quency, and plot starting frequency are printed in the upper right-hand 
corner of the plot. An example of power spectrum plot is presented in 
Section 3.4. 

3.2.3 Tape Output 

SAP generates the SAP plot tape and the SAP spectrum tape. The SAP 
plot tape contains one data file consisting of two records. Record #1 
contains 1) the upper/lower plotting surface designation. 2) case title. 

3) plot frequency spacing, 4) carrier frequency, 5) plot starting frequency, 

6) number of spectral points to be plotted, and 7) a SAP internally generated 
parameter (IFC) related to a centered spectrum plot. Record #2 contains up 
to 480 spectral values (in decibels). The SAP spectrum tape contains one 
data file consisting of the signal processed through entry point 10. This 
file contains NCIDT record pairs with NRIDT words per record. A record pair 
consists of two records - the first containing the real part of the signal 
and the second containing the imaginary part. 

The above two tapes are generated by a normally terminated SAP execu- 
tion. At intermediate stages of signal processing, the output from a given 
module is stored on tape as a single file consisting of NCIDT record pairs 
and NRIDT words per record. As each entry point -is passed during execution, 
a message is printed indicating the number of columns of output data (NCIDT) 
and the number of rows of output data (NRIDT). 
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3.3 OPERATING PROCEDURE 


The operating requirements and procedure for executing SAP is presented 
below. 


I. REQUIREMENTS 

A. Program Deck 

B. Input signal data tape mounted on tape unit IDT. If signal 
is to be generated by TSGEN a spare tape is mounted on tape 
unit IDT. 

C. A spare tape is mounted on tape unit 1ST. 

D. CS-1 subroutines MAGTAP and CL0CK on paper tape. 

II. OPERATING INSTRUCTIONS 

A. Load program deck 

B. Set sense switch 6 up 

C. Compile program 

D. Check load map for starting core storage locations of 
MAGTAP and CL0CK 

E. Load CS-1 routine MAGTAP (paper tape) into correct 
core location 

F. Load CS-1 routine CL0CK (paper tape) into correct core 
1 ocati on 

G. Reset sense switch 6 

H. Execute program 

III. OUTPUTS 

A. Computer printout 

B. SAP plot tape 

C. SAP spectrum tape 

IV. RESTARTS 

A. Check printout for restart information. This information 
will include: 

(1) Last successful entry point completed 

(2) Data tape unit 

(3) Spare tape unit 

(4) Number of rows of data on data tape 

(5) Number of columns of data on data tape 
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B. Reload data cards with the above changes, also load ENDATA, 
SYSEND, and blank cards. 

C. Continue operation starting with Operating Instruction 1 1 - H . 

V. ADDITIONAL RUNS WITHOUT RECOMPILING 

Additional runs may be made by 1) loading data cards, SYSEND and 
ENDATA control cards, and 2) starting with Operating Instruction 
II-H. 

3.4 SAMPLE CASE 

A sample case is presented to demonstrate the program capabilities 
and options. The sample case includes specification of the problem, sample 
coding form, sample listing of. the program input and output, and the associ- 
ated power spectrum plot. 

Example . Compute the power spectrum of a 1 MHz carrier (FC = 1 x 10^) 
frequency modulated (IFM = 1) by a square wave of frequency 1.5625 Hz 
(TP = 0.639990) and amplitude unity (A = 1.0). Exactly one period of the 
square wave will be generated internally (ISOPT =3). A peak frequency 
deviation of 200 Hz (BETA = 200.0) is desired. No filtering is to be per- 
formed. The output power spectrum, centered about the carrier (FSTART = 

-1.0), is to be both printed and plotted using the maximum available fre- 
quency resolution (NAVE =1). The spectral region of interest is 
1 x 10 6 + 375 Hz. Based on analytical considerations, the amount of 
aliasing is estimated to be small in the region of interest if the number 
of samples is selected to be 1024 (NRIDT = 32 and NCIDT = 32) and the 
sampling time interval is selected to be 0.625 x 10" 3 seconds (DT = 0.000625). 
(The validity of this estimation may be experimentally verified by comparing 
with a second power spectrum computation in which the number of samples is 
doubled and the sampling interval is halved. The estimate is valid when 
the agreement between the first and second computations is within an 
acceptable tolerance in the region of interest.) 

3.5 TIMING DATA 

Typical timing data for SAP is presented in Table 3-4. 
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♦**♦•***•***•***•***»»*♦»******••******»****-•*«-*»*«» 


IDT I£P IF,/. IFOPT IFTYPI IFTYPZ" IPLPOS" IPMoOE" 150PT - r5T ”KSK _ IiCIDT NCI NP1 NfZ NRIDfWrNAVE IZPRT 

I * I 0 0 0 0 0 3 20 32 000 32 0 1 000000 

b£TA = 0.2OuUUO00E03>“ ~ ' ' 


A ^ rtA ” Oi_uo> bETA = 0.2OuUUO00E03» OELT = 0.6.Y500000E-03. FSTART = OEOO. FC = 0. 10000000E07 

***** ** *** ********** ********************************************* '/♦♦***********#********#******#**************************ti******** 

TP A 

0 .63999990EOQ 1 . OOOOOOOOEOO 



“** PRogkAm has passed entrt point Tool m 


** DATA IS_CURRENTLY ON UNIT 1. __ ** 

** spake Tape is cmmENtLY on unit z '.' **' 

«« Ther e AHE 32 ROWS AN O 32 COLUMNS OF DAT A. »* 
*» EXECUTION time = 0 MINUTES* 2.'968 SECONDS »* 

♦ * **. 

***•*•*••**•*•****•**•*•'>••*•••***•***•*••••••**•• 




** *« 

** "PROGRAM HAS PASSED ENTRY POINT 200. ** 

“ DATA IS CURRENTLY ON UNI T 1. ** 

** spare tape is currently on unit 2. •* 

** there are 32 rows and sz columns of data. »» 

*« EXECUTION TIME = 6 MINUTES* 4.538 SECONDS ** 

«» »• 


******************* ********************** ******** ** 


♦♦•♦♦•♦•♦♦•♦••••♦••♦♦♦•••I************************* 

*• 

** PROGRAM HAS PASSED ENTRY POINT 300. *• 

«« DA T A IS CURRENTLY ON UNIT 1 . ♦ * 

♦* spake tape is currently- on unit z . »* 

»« THERE AKS 32 ROWS AT ID 3 2 COLUMNS OF DATA. »♦ 

** EXECUTION TIME = 0 MINUTES* 16.516 SECONDS ** 

.** , * * 

*********** ***********»*“*******i****4************ 



************************ «************************^ 

♦* 

** PROGRAM HAS PASSED ENTRY POINT 400, *• 

*♦ DATA IS CURRENTLY ON UNIT 2. ** 

** spare tape is currently on unit i. ** 

** there are 32 rows and 32 columns of d at a. « « 

*• EXECUTION TIME = 1 MINUTES. 3.427 SECONDS ** 

<• *• 


♦ ♦♦A********************.'************************* 


********** **•**•**•***•**••* ************* ********* 


*♦ _ ** 

** PROGRAM HAS PASSED EnIrT POINT '500\~ •* 

** DATA I S CURRENTLY ON UNIT 1. •< 

** SPAKE TAPE IS CUFIRENTLT on UNIT 2. ** 


• * there are_ 32 rows_ and 32 columns of data. ** 
*« EXECUTION TIME = 0 MINUTES. 9.347 seconds ** 

•* _ ** 

*********************************** *.»*****««**»** 


3 

3 

5 
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******************* *********************** ******** 
** **_ 

_ - pKOGhAM HAS PASSED ENTRY POINT 600. ** 

** DATA IS CURRENTLY ON UNIT 2.__ _ _._** 

** SPARE TAPE IS" CURRENTLY ON UNlt 1> ** 

** THERE ARE 32 ROWS AND 32 COLUMNS OF DATA..** 
** EXECUTION TIME = I MINUTES* 1.521 SECONDS ** 

** ** 

»*-********************* *************************** 


************************************************** 

** _ ** 

** PROGRAM HAS PASSED ENTRY POINT 700. ** 

** DATA IS CURRENTLY ON UNI T 1. i* 

** SHAKE tAPE IS CURRENTLY ON UNIT 2. ** 

*♦ there are 32 rows and 32 columns, of data. **. 

** execution time ="o minutes* 26.222 seconds ** 

** ** ... 

“■ — ************************************************** 

— ** ************************************************ 

** _.** ... 

** PROGRAM HAS PASSED ENTRY POINT 800. ** 

** DATA IS. CURRENTLY ON UNIT_2. ** 

“■'** ~~~spare Tape is currently on unit 1 . ** 

** THERE ARE 32 ROWS AND 32 COLU MNS OF DATA. ** 

** EXECUTION TIME = 1 MINUTES. 8.054 SECONDS ♦* 

* * _ * * ... 

" “■***'********************■ 1 ************************** 


************************************************** , 

*♦ ' ** 

— **~ PROGRAM HAS PASSED ENTRY POINT 900. ** 

*♦ DA TA I S CURRENTLY ON UNIT_2. .**.._____ 

** SPARE TAPE IS CURRENTLY ON UNIT 1. ** 

«» th ere ar e 32 rows a nd 32 columns., of. data •_** 

** EXECUTION TIME = 0 MINUTES. 0.616 SECONDS ** 

>« ** 

************************************************** 


************************************************** 


*♦ _** 

** PROGRAM HAS PASSED ENTRY POINT 1000. ** 

** DATA IS CURRENTLY ON UNIT 1. f*. 

♦* spare Tape is - currently on unit 2 . ** 

** THERE AHE 32 ROWS AND _32 COLUMNS OF DATA. ** 
- ♦* EXECUTION TIME = 0 MINUTES. 6.806 SECONDS ** 

** ** 


— — — — — • — ************************************************** 

*♦««*♦♦»»»*»**»♦»«♦ ♦♦*♦♦♦*♦***«♦*«♦»♦»*««*♦««♦*»«♦♦*»*♦***♦♦♦***»************♦***♦* ********************************** ^-t*«»******«*A 


PRINTED IN U S A. 
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u.999703l2EUb 
u.9997U4o9tbb 
j;9S97u6250b 
u . 9997u7dlb Ub 
' u . 99970937608 
d • 99971094E0O 


7u00uUOUEU2 

43B09191E02 

7 b 0 u 600 d £02 

4354o8 79E02_ 

7lu6uu6uEu2 

4327abUlE02 

7oPooOOOE02 

43003344E02 

/0000ddoE02 

42 7214o4Eo2 

7u0060duEU2 

4243198oEU2 

7oOOUGOO£o2~ 

42134b45E02_ 

7U0004UUEo2 

41d2b7o7EU2__ 

7bU0U0UUtu2 

4lblo28b£u2 


0 . 99985312E06 
0.99985469606, 
U • 999b5b25E06 
0^999u578lEqb 
' 0.99985937608 ' 
0.999o6U94E06 


• 70000000E02 
i_32229831£02 
.70000000602 
■ 32738255602 

• 700Q0O0OE02 
. 33210198E02 


0.10000000E07 
0.10000000E07 
U.10000000E07 
_o.ioogopooE07 
O.IOOOOOOOEO?' 
0.1 OO OOQOOFO I_ 


70000000E02 

40261462E02_ 

70000o60Eb2 

4023949QE02 

70 0000 0 0E02 " 

4 0213045602 


0.10001937E07 -0.70000000E02 

‘ 0.1000 195 3 607 -0. 13 374B 3 3F02 


0.10001969E07 -0.70000000E02 

0.10001 964EQ7 -0 .38930 1B3E 01 


0.10002000E07 
0.1 000201 6E 07 


'6799W6250E06 
0.99986406608 
0.99986562EO6 
0 .999d6719E06_ 
b . 999d6875E06 


-0.70000000E02 

-0.33648802E02 

-0.70000000102 ' 

-0.34 0583S4EQ2 

-o.7bbobbb"oE02 

-0.34441582E02 


0.10000312E07 
_0. 1000.0 328E 07 _ 
0.idobb344E07 
_Q • 100Q0359EQ? 
0.10000375E07 
0 . 10000391F07 


0.10000406E07 

0.10000422E07_ 

o.iobob437E07 

0.10000453607 


70000000E02 

401B2116E02 

70006000E02 

40146590E02 

70oob6ooE02 

401 06512 F02 


0E00 

=fiii9529M3£lll_ 

0.10002031E07 — 0.70000000E02 

0.10002047E07 ^0,1555474215? 

6.lbbb2062E07 -0.70000000E02 

0.10Q02076E07 -0. 1S05Q403E02 


0.10002094E07 -0 
0. 10002109F0 7 -0 
0.10002125E07 -0 
0.10002141 E07 -0 


-0.7000000 0602 
-0^34801197E02_ 
-0.70000000602 
-0.35139125602 
-0.70000000602 
0 . 999 d 7 0 3 1 E 0 6 -0.35457416602 


0.10000489E07 

0.10000484E07 


0.99987187E06 

0.99987344E06 

0.99967500E06' 

U.999d7656£06_ 

0.99987812606 

0.99987969E06 

0.99986125E06 

0.999d828lE06_ 

0 .99988437£b6 

0.99988594E06 

0i99988750E06 

0. 9998890 6E06 


-0.70000000602 
-0 • 35757532602 
-0.70000000602 
-0.36041 2 13602_ 

-o.7bobbbboi:o2 

-0.36309248602 


0.10000500E07 
0. 10Q00516E07 
0.10000531E07"" 
0.10P00547E07_ 
0 ■ 10000562E07 
0. 1 000 057BE 07 


70000000602 

40061752EQJL 

70000000E02 

40012291E52_ 

70000000E02 

39957974F02 


0.10002156E07 -0 


7U00uUU0EU2 

411923U4E02 

70OOUOO0EO2 

40859211602 

/u0uUUUUt02 ’ 

4051b54otU2 


u.94971«.bOLUo 

0 . 999714UULU6 

' J. 99971bd2Lu6 
_u.99971719LUd 
u.9997lo7bEuo 


7UOOUOOUEU2 
4UldU9b4EU2 
7bOOUOOuEU2 ” 
39794130602 
70bUuU0UEU2 


0.99989062606 

0.999d9219E06 

0 . 999d9375E06~ 

0.999d953lE06 

0 . 99989687EU6 " 

0.999d9844E06 

0.99990000E06 

0 . 999901 566 U6_ 

b. 99990312608 ' 

0.99990469E06_ 

b.99990625E”b6 


=0770000000602 

-0.36563057602 

-0.700000066:02' 

-0.36803465602 
=0.70000000602 " 
-0.37031491602 


0. 10000594E07 
0.10000609607 
0. lbobo625E07" 
o,igogg64iEP7_ 

' 0 . 10000656E07 
0 . 10000672E07 


70000000E02 
398988 13£0?_ 
7b0OOOO0E02 
39B34654E02_ 
70000000E02 
jg765a.I7j;02 , 


0. 10002187E07 -0 

-0.100022D3EQ1 -0 

0 . 1 00022 19E07 -0 

. 0 , 

0.10002250607 -0 

_a,10P02266EQ.7 . . -Q 


70000000E02 

2103114QEQ2 ' 

.70000000E02 

2J27UI13ED2 

70000000602 


-0.70000000602 

-0.37247715E02 

=b.7oobboooEo2 

-0.37452973E02 


0 . 10000687E07 
0.10000703607^ 
0. 10000719E07 
0. 10000734E07 
'b.i0000750E07 
0.10000766E07 


70000000E02 

3969O9J2E0g_ 

70000000E02 

3?611222E02_ 

7b66ooooE02 

3952 6010E02 


0.10002281E07 -0 

_g.l0002297EQ7_-fl 
0.10002312E07 -0 

0.1000232BEQ 7 -0 

0.10002344E07 -0 


70000000E02 
. 265 7866 3£02_ 
70000000E02 
27077466602- 
70000000E02 
29M9it63EH2- 


70000000E02 

39435226E02_ 

70000000E02 

39338624£02_ 

70000000E02 

39236144F02 


0.10002375E07 
_g.lO0U2391E0Z_ 
0.10002406E07 
- 0,1 
0 . 10002437E07 
. P ■ 1 0P529.53J60.7_ 


70000000E02 

*30040235602- 

70000000E0? 

30?&J4BDED2_ 

700O0000E02 

■3ian7l?5En2 


0.10002469E07 -0 

-£L,i 0502484607 -0 

0.10002500E07 -0 

0 .100 02bl6E07 -5 


3.70000000E02 
* 325fl4Q71E0iL 
•70000000E02 
533049ol3£02_ 
.70000000E02 


-0.70000000E02 
-0 . 376476 13E 02 


-0.70000000E02 
-0.37833180E02 
=o.7booboboEb2 ^ 
-0.38008698602 
-0.70000000602 


0.10000781E07 
0. 10000797E07_ 
' 0.i00008i2E07 
0.10000828E07- 
" 0. 10060844E07 
0.100008 59 E07 


70000000E02 

39127507602 

70000000E02 

39012593E02_ 

70000000602 

3889Q980E02 


0.10002531E07 -0 

-0,108 02 54 7E8I -0 


0. 10000875E07 
0 . 10000S91E07_ 
0.10000906E07 
0.10000922E07 
0.1C000937E07' 


70000000E02 

38762703602 

7boobbooE02 

38627293602 

7000bobbE02 

38484556 E02 


0.10002562607 -0 

0.10002578EQ7 _-g 
0. 10002594E07 -0 

. 0 • 10002609E07 -0 

0.10002625607 -0 

0 . 1 0002641F07 -0 


70000000E02 

39608.035602 

70000000602 

,35201571602 

.70000000E02 


70000000E02 

38334061602 

70000000602 

38175607E02 

70000000602 


0 • 1 000265SE07 -0 

0,1000267260? -0 

0 • I 0002687E07 -0 

0.10002703607 rO 
0. 1000271 9E 07- -0 
J>, 10002734F07 


•70000000E02 

,36294761602- 

•70000000E02 

.36801145602- 

.70000000E02 

372fi3?.Q2EQ2_ 


0.10002750E07 -n 

. 0.10002766607- -0 

0.10002781607 -0 

. 0 • 10002797607 -5 

0. 10002812E07 -0 


.70000000602 

37745879602 

70000000602 

38188007602 

70000000602 


.70000000602 
. 39020959fj)2_. 
.70000000E02 
39414584E02- 
.70000000602 
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«i Id7t.uo 

J. WS7<:j‘4he.uo 

-wiV4l4bb3Eu2 

0.99990781E06 

-0.38175606E02 

0 • X0000953E07 

-0 . 38008700E02 

0.10002828E07 

-0.39794136E02 

-u. ?UU0UUDUEtl2 
-U.3902 U^dOEU2 

0.99990937E06 

0.99991094E06 

-0.70000000E02 

-O.30334D61EO2 

0.10000969E07 

0.10000984E07 

-0* 7000000 OE 02 
-O.37833180EO2 

0.10002844E07 

0.10002859E07 

-0.70000000E02 
— 0.40160956E0? 


j.99972o5oloo 

u.9yy7^bliEtiii 

o.9997<:9o9c.0b 

O.9997oi25c.oo 

u.yyy7j£oit.uo 

0.9997j4j7tLUo 

0.9497bb94Luo 

ui9997i7i)Ucuu 

0.9997j9docou 

~o.99974Uo2c.ub 

0.999/4<.i9c0b 


“0 • 3obl2857Lu2 
-o. /uouuuuuEoa 

“!>• Jaloo0UoE 02 
-U. 7uuuuuuuEU«f 
-0.37745o78E02 


”379997^375^00“ 
J • 999745 JIl.06 
J*9997hoo7£u6 ' 
uu9997-*d44cbb 

~u .999 76uuOcuo ' 
0.9997b! but 00 


-0. 7U0UUUU0EU2 
-0. 0728o902E02 
-0. ?0003uuutu2 " 
-0. Jo80il44E02 
7uOuuU0UEO2 
’0 . 3o29h 70tE02 


0.999754o9tuo 
~u799975625EGb 
u. 999707a it. 00^ 
0.999709 J7tuO 
0.9997oo94cub 


~=u~. 700033 
- 0 . 357otbbUE02 
-0. 7uOuuOuu£02 
-0. J5201074E02 
■-0.70UdU00uEd2' 
-0.34608031E02 


0.99991406E06 
0.99991562E06 
0.99991719E06 
b .99991 87 5£ 06 
0 . 9 9992031 E06 
'0 
0. 
o 
o 
b 
o 


99992187E06 - 

99992344E06 


-0.3 S 4a4S S7E02 

-0.700U0OO0E02 

-0_.38627293E02_ 

-U.70000000E02 

-0.38762703E02 


999920bOEb6 
9999b656E06 
99992312E06 -I 
999929b9E06 


0.70000000E02 
(K 3889098 1E02_ 
0.700000b0E02 
0 ■ 390 12591E02 

o.7ooooaooE02' 

0.39127S07E02 


0.10001000E07 
0. 10 00101 6E07_ 
'b;ioooib3iEb7 
0.10001047E07. 
b.i0001062E07 
0. 10001078E07 


“0799975753E3Sr 
u.9997o4UutU6 
U . 9997o562cG"6~ 
0.9997o7 1 9EUO 

0.9997ob7bt0fa 

O.99977031tOb 


-0 • 7uOOu030E02 
-0.339772O4E02 
-0. 7u0duO0OEu2 
-U.33304901EU2 

-O77oouuuodtu2 

-0.32584U71E02 


99993281E06 
99993437E06 ”-i 
99993094E06 
999937S0E05 
999939U6E06 - 


U.99977344cUb 
~3.99977533Eu6 
U.99977obbLUo 
”o.999 77b I'EtTob '' 
0.999779o9E0b 


-O.7uOOu3o0E02” 

-0.31807123E02 

'=u.7uduuOdoEu2” 

-0. J09bj4UUEU2 
-0. 7UOUU00OEO2 
■U.3004U2 odE02 


r99994362E06 

. 99994219E06 

.99994375E06T” 

.99994S31E06 

,99994o87E06 

.99994844E06 


fl.?booOo“ooEoi 
0.39236144E02 
0 -70 0 0 0 0 0 0E0 2~~ 
0 . 39338 625EQ2 
0.70000000E02 
0.39435225E02 


0 . 1 0 00 1 094E07 

_0 . 10001 109E07 

0.100011 25E07 
J).10001141E07_ 
0 . 1000 1 156E07 
0. 10001172E07 


-0.70000000E02 
- 0.37647 6 0 4E 02 
— 0.70000000E02 
j-0.37452971E.02_ 
-b.70600660E02 
-Q. 372477 14F02 


-0. 7000000UE02 
— 0.290 19464E 02 
-0. 7uOOOOOUEfl2 
-0.27B77464E02 
-0. 7U00U0U0EU2 
-0.2b57bo62E02 
-0. 70000000E02 
-0.2b071 J82EU2 
-0. 7oOOUOOOEU2 
^O.c3271413Eu2__ 
-0.' 70Q0OOUOE02 
-U.21031139E02 


0 . 7b00b000Eb2 

0.39526011E02 

b.7obbbobOEo2‘ 

0.39611223E02_ 

0.70000060E02 

b.39690916E02 


b. 10001187E07 
0. 10001203E07 
0.ib001219E07 
O i 10001234E07_ 
~0.1000i250E07 
0. 10001266E07 


-0.70000000Eb2 

-0.37031491E02_ 

-b.700ddbbbE02 

-0.36803463E02 

-0.70000000602 

-b . 36563056E02 


0.10002875E07 
0.10002 891E 07 
0 • 1 OOb2906E07 
0.10002922E07 
O.Io002937E07 
Q.lb002953 Eb7 


^9953doEo6 -d.7ao00000£02 
99995156E06 -IK 39765474E02 
> 999953 I2£fl6”~ — 0 .7bObObOdE02 
9999S469EU6 -0 . 39834654E02 
99995625E06 -0.7bbQbObOEb2 
99995781E 06 -0.398988 13E02 


0.10001281E07 
_0. 10001297E07 
0.10001 3 12E07 
0,1000132 BE 07_ 
0.1066l344E07 
0. 10001359E07 


-0.70 00000 OE 02 
— 0_» 36309248£02_ 
-0. 7000000 OE 02 
jjO . 36041214602 
-O.7O060b6bE02 
-0.35757 5 30E02 


0.10002969E07 

0,10002964E07_ 

d.io003000E07 

0. 1000301 6E07_ 

0.10003031E07 

0.1 0 003Q4 7E07 


-0.70000000E02 
-0.4 0515 546E02 
-0.70000000E02 
-0 .40859 21 2E02 
-0.70000000E02 
-0.<H192?06ED2- 


W99593 7T06 -0.700000 0 OE 0'2 

99996094E 06 -0 . 3 9957974E02 

99996250E06 = OT;70000CObE02 ' 

99996406E06 -0.40012281E02 


0 . 1000 1375E07 
_0.10001391E07_ 
b.idbbi406E07 
0.10001422E07 
0.100oi437E07 
0. 10001453E07 


•0.70000000E02 
•0.35457415E02 
■U.70000000E02 
■0.35139125E02 
■d.7oooddooE02 
— 0 . 3480 1 197F02 


0.10003062E07 
0. 10003078E07 
0.10003094607 
0 ■ 10003109E07 
0.1000M25E07 
0 . 1 0003 1 41E07 


-0.70000000E02 

-Q^S1516294ED2_ 

-0.70000000E02 

_rO i 41#28780E02_ 

-0.70000000E02 

J^.42i34844Ea2_ 


99996562606 -0 . 70000C00602 

99996719E06 -0.40061753E 0 2 


0 •999968?5Eb'S” 
0.99997031E06 
0.99997r87E06 " 
0.99997344E06 
0 . 99997500 £06 
0.99997656E06 


0. 10001469E07 
0 ■ 1 000 1484E07_ 
b.idodi5odEd7 
0jl00_Q1516E07_ 
'o.ioodi53iE07 
0. 10Q01547E07 


— 0.70000000E02 
-0^3444 158 IE 02_ 
-o.7dooooodEd2 
j?0_.34058354E02_ 
-d.7dooooodE02 
-0 , 3364880 2E 02 


0 • 10003156EQ7 
0.10003172E07 
0.10003187E07 
_0.10003203E07_ 
0 . 100 0 32 19E07 
0 . 10Q03234E07 


■0 . 70000000E02 
— 0.424 J1984EQ2_ 
-0. 7000000 OE 02 
r0.42721480E02_ 
-o.7odoboodE02 
-0.43003343E0 2 


-0.70000000E02 

-0.33210202E02_ 

-d.7doboooOE02 

-0.3 2738251E02 


0. 10003250E07 
Q.1 0003266F07 
0 » 100 0320 1E07 

_0_.10003297E07_ 

0.100d33i2E07 

0.10003328E07 


0 • 7000 0000602 

•0j.4327 8 5aiEa2_ 

•0.70000000E02 

r0.43546860EO2_ 

-0.70000000E02 


— 0.70000000E02 
-0.32229830E02 


0 . 1 0003344E07 

0.10003359E07_ 

0.10003375E07 

0.1000y91E07_ 

d.lddd3it06E07 

0.10003422E07 


-0. 7000000 OE 02 

J-0..4H06H893EQ2 

-0.70000000E02 

^Oj44315454F02 

— 0.70000000E02 


-0.70000000E02 

j:0..‘l4e00925EB2 

— 0.70000000E02 
-0.45034448E02 


-0.70000000E02 


J.9997ai25E06 

0.9997a^olL0o 
”OT9997o437£o 6 
0.999 7ob94t. Ob 
o.999787but.0b" 
0.9997bVuoEUo 


-0.70000000E02 
-O. 401 0651 3E 02 
-o.7oOoooooe62' 

-0.40146592E02 


J . 999 79062Eoo 
0.99979219E06 
-•3799979375E05 
0.97979531EU6 

””3r999796ar£o6' 

0.9997 9O49EUO 

0.9990olbot.oo 

0i999803l2t0o 

0.999d04o9EU6 

0.9998O625E06 

0.999807alt0o 


-0. /OOOOU30E02” 
-0.18Q50402EU2 
-b. 70003000102 
— Q.135547 h1Eu2 
-u;7o0000ooE(j 2” 
-0.3952960dE01 


0.99997969E06 • 

”3. 999981 2 5£ 06 — ■ 
0 . 9999828 1 E06 
0.3999m 43 7E06T 


-0.70000000E02 
-0 . 4Q1821 15E02 


■3TT4'2'6370bE-06 
-Q.3o9301db£01 
-0. 7UOOOOOOE02 
-0. 13374833EU2 
-0. 7u000000t02 ' 
-0.17750607E02 


0 . 7oOoOOOoE02 07 

0^402 13046E02 0. 

3.7bOOOCOOE02 0. 

0_.40J39491E02 0 . 

0.70000000E02 0. 

0.99998594E06 -0 . 40261464E02 

”07 

U 


10001562E07 

10001578E07_ 

10001594E07 

10001609E07_ 

1000 1 625E07 

10001641E07 


0.9998u937tub - 0 . ?000o00UE02 

_U.9990lU99k.Ub — 0.2061121 5E 02 

*************** ******************** 
*********************************** 

Kb A* =‘ ” 0T2bu33932E0U 


”0';9999 fl 75oE06 - 

0.99998906E06 - 

0 . 99999062E06 - 

0. 999992 19E 06 

O. 99999375 E 06 =i 

0.99999531E06 


37700000 00 E02” 
0.40279015E02 


1000 1656E07 

10001672E07_ 

10001687E07 

10001703E07_ 

lbobi719E07 

10001734E07 


— 0.700 0000 0E02 
—0 .31679006E02 
-07706d0000E02 
-0.31079678F02 
-o.7obodoooE02 

-0. 3042 366 7E 02 


0.999996B7EOb - 

0. 9999964 4E 06 - 

**************** 

**************** 


T77oooo6ooe62 0. 

0^402921 l 4tiE02 0. 

o. 
o. 


0.70000000E02 
0^40300 a 97E0 2 


0.70000000E02 u 
0.403Q5273E02 • 0 

♦******♦•*•* ******* 
******************* 


1 000 1750E07 

10001766E07_ 

lbdbiteiEof 

10001797E07_ 

i00018i2E07 

10001 828E07 


-0.70000000E02 

-0.29701113E02 

-0.70000000E02 

-0.28898589E02_ 

■0.70Odo"000E02 

•0.279986B0E02 


0.1000 34 37E07 
0. 10003453Fn7 
b.iodb3469E07 
0. 10003484E07 
0.10003b0OE07 
0.1 00035 1 6Ej} 7_ 


10001844E07 

10001859E.07 

************* 

************* 


-0.70000000E02 
_rL*^ 97 6662E02 
-0.70000000E02" 
^0^25798450602 
-0.7O0d6b0OEO2 
-0.24 4 11345E0 2 


0 . 10003531E07 
(K 10003547E07 
0.10003562E07 
0.1000357BE07_ 
b.i0003594E07 
0 . 1 0003609E07 


-0. 70 00 000 OE 02 
^,95489462E02_ 
-0.70000000E02 
-0.45 71 062 2E02 
— 0.700 0000 0E02 
-0.4 592617 4E02 


•0.70000000E02 
r0.46l38927EQ2_ 
-0.70000000E 02 
-0.46347262E02 


-0.70000000E02 
-0_. 22731 645E02 

**•*«••*•»••**• 


0 . 10003625E07 
_0.10003641E0.7_ 
0. i0O03656E07 
0.10003672E07 
0. 10003687E07 
0. 10003703E07 


-0 • 70000000E02 
-0 .46551933EQ2 


-0.70000000E02 

^P.*46752597E02_ 

-0.70000000E02 

70 s 46959Pt4£S2_ 


0.10003719E07 
Oj 10003734107 

****••****•*••*•• 

***************** 


-0.70000000E02 

rg, 47 14 J t8X5EQ2.. 


-0.70000000E02 

rO,_HI3J4567£02- 
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*************************************************** ********************************************************* **+**?±* ****** *****•***. 

' rtUN Tt^MINATfcD. PLOt DATA TAPE' GEfJERATED ON UNIT 2 

*********************************************************************************************************************************** 

_ EX£C0TI0N TIME = r MINUTES. 37.985 SECONDS ** 

** ** 

~ ♦*#*******♦************;****♦***********♦♦********♦ 







TABLE 3-4 

TIMING DATA FOR UNIVAC 1230 SPECTRAL ANALYSIS PROGRAM 


NUMBER OF 
SAMPLES 

1,024 

4,096 

4,096 

16,384 


NCIDT X NRIDT 
32 x 32 
64 x 64 
64 x 64 


EXECUTION 
' TIME 

5 min, 38 sec 

20 min, 19 sec 

32 min, 13 sec 

80 min, 41 sec 


COMMENTS 

No filtering performed 
No filtering performed 
Bandpass filtered 
No filtering performed 


128 x 128 



4. PLOT GENERATION PROGRAM 


4.1 DESCRIPTION | 

The PLTGEN Program is a stand alone FORTRAN program, designed to read 
a data tape generated by the SAP program, convert the data into plot commands 
for the EAI 3440 DATAPLOTTER, and write these commands on a magnetic tape. 

The plot subroutines may be divided into two categories: the basic 

plot routines and the routines designed specifically for the construction 
of the spectral analysis plot. 

The basic plot routines provide the general functions of a plot package. 
These subroutines initialize the scale factors, raise and lower the pen, 
move the pen from one position to another, write formatted numbers, print 
character strings, and write these formatted plot commands on a magnetic 
tape. 

The special plot routines utilize the basic routines in forming the 
grid, annotating the axes, and plotting the spectrum. The main program 
PLTGEN, is the driver, which calls the special subroutines in the proper 
sequence to complete the desired plot. There are nine subroutines in addi- 
tion to the main program. The basic subroutines are: COMAND, INITIAL, 

NUMBER, PEN, PLOTPT, and PRINT. The specific routines are ANNOTE, PLOTER, 
and GRID. These subroutines are each briefly described below. 

COMAND - Subroutine COMAND is called each time a plot command is written 
on the output plot tape. An X and Y command code, which specifies the 
comnand function, and corresponding X and Y numerical value is received 
through the calling arguments, formatted into an eight word record, and 
written on the output plot tape. 

INITAL - This routine is used to initialize the EAI 3440 Dataplotter by 
constructing the scale factor, data offset, and board offset comnands. 

These three commands must be the first three written on the output plot 
tape. 


4-1 



NUMBER - Subroutine NUMBER converts a floating point number into a formatted 
Hollerith string. The Hollerith string is contained in an array having one 
Hollerith character stored left adjusted in each word. The two formats 
available are the F7.1 and the scientific notation E10.3. The formatted 
string is then passed to subroutine PRINT for further processing. 

PEN - This subroutine is used to select the plot pen or the printer. After 
the plot pen has been selected, subroutine PEN is used to place the pen in 
the up or down position. 

PLOTPT - The purpose of subroutine PLOTPT is to construct a command to move 
the plot pen from its present position to a new (X,Y) position. The plot 
pen may be in either the up or down position when the PLOTPT command is 
given. 

PRINT - Subroutine PRINT converts a string of Hollerith characters into the 
EAI 3440 DATAPLOTTER character set. The Hollerith string of characters 
must be contained in an array having one Hollerith character left adjusted 
in each word. The Hollerith string will be printed with the first character 
at a specified (X,Y) position. The orientation of the characters and the 
distance between characters must also be specified through the calling 
arguments. 

ANNOTE - Subroutine ANNOTE calculates a scale appropriate for the amplitude 
and frequency bandwidth of the spectrum being plotted. The carrier frequency 
is used as a bias frequency and each frequency label is biased by that 
amount. The placement of the tic marks along the x axis is determined so 
that the corresponding labels are in even increments of frequency. The Y 
axis is labeled from 0.0 to -70.0 dB using increments of -10.0 dB. A plot 
title and other pertinent spectral frequency information is printed in a 
title box in the upper right corner of the plot. 

PLOTER - This subroutine is used to draw' the spectral lines on the plot. 

The input data array is received through Labeled Common /A/. A straight 
line is drawn from the x axis to each amplitude value in the data array. 
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GRID - The purpose of this routine is to draw a rectangular box 25 inches 
by 7 inches with tic marks 1 inch apart along the two 7 inch sides. A 
small title box is drawn in the upper right corner of the rectangular box. 

4.2 OPERATING INSTRUCTIONS 

The PLTGEN Program can be loaded into the UNIVAC 1230 computer by 
reading the card deck or the magnetic tape containing the compiled program. 
Execution of the program requires a small deck consisting of five cards 
(see Figure 4-1). These cards must be read for each data case. These 
cards will be described in more detail below. The PLTGEN Program requires 
approximately five minutes to load from the card deck and about one minute 
to read from the compiled tape. The execution time required for one data 
case is approximately one minute. 

In executing the program, the UNIVAC 1230 operator should be instructed 
to load the PLTGEN Program from the desired input source. He should also be 
instructed to mount the SAP output spectrum data tape on Unit 2 and a blank 
plot tape on Unit 1. The five cards required for executing the program 
should be available for the operator at this time. Multiple cases may be 
stored on one plot tape by mounting different spectrum data tapes on Unit 2 
and re-executing the program. The PLTGEN output plot tape on Unit 1 should 
not be rewound until all desired cases have been run. The plot tape may 
be plotted at this time or stored by the operator for plotting at a later 
time. 

Before plotting the PLTGEN output tape, the operator should be instructed 
to make the proper switch settings on the EAI 3440 plotter control panel. 

The switch functions and proper settings on this panel are discussed on 
the following pages. If multiple cases have been written on the tape, the 
operator should be informed how many plot files are on the tape. 
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Figure 4-1. Execution Card Deck 
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4.3 OPERATING CONTROLS AND INDICATORS 


a. Control Panel 20.700 

CONTROL FUNCTION 

X AND Y CHARACTER 
SELECT 


REQUIRED 
SWITCH POSITION 


SIGN This rotary selector switch allows the 11 

operator to select any of the twelve 
characters in a word as the sign charac- 
ter. With this switch in the zero posi- 
tion the system reads all input data as 
positive. 


COMMAND 


THOUSANDS, 
HUNDREDS, TENS, 
AND UNITS 


The setting of this switch determines 1 

which character in the word functions 
as the X or Y command. If either switch 
(X or Y) is rotated to zero, both con- 
trols are considered to be zero, the 
command codes from the tape are ignored, 
and all data is plotted. 


These 13-position selector switches 
provide the operator with a means of 
selecting any of the twelve characters 
in a word as the X or Y thousands, 
hundreds, tens, and units value. 


Thousands 2 
Hundreds 3 
Tens 4 
Units 5 


Y 

2 

3 

4 

5 


X AND Y WORD SELECT y 

TENS AND UNITS The X and Y word selectors allow the Tens 0 0 

operator to select any of the first Units 1 6 

159 words in a record as the X and Y 

word for processing. The TENS switch 

functions as a XI 0 multiplier and the 

units selector setting adds to the TENS 

switch setting. These switches select 

the word according to the setting of the 

WORD SIZE switch. Therefore, if the 

words selectors are set to 150 and the 

WORD SIZE switch is on 10, the first 

character in the selected word would 

be the 1 49s t character in the record. 


4-5 



CONTROL 
WORD SIZE 


TAPE MODE 


PLOT MODE 


MAX POINT 
DISPLACEMENT 


FUNCTION 


REQUIRED 
SWITCH POSITION 


This 11 -position rotary selector is 12 

used to set the number of characters 
in a data word. The word size is 
variable from 2 to 12 characters. 

This switch determines the speed at which Auto 
information from the tape is read and 
plotted and also controls the tape trans- 
port. In the AUTO position, the tape is 
momentarily stopped at the end of each 
record and the command codes cause the 
appropriate action to occur. In FREE RUN, 
the tape moves at 75 inches per second and 
the tape does not stop unless a command 
other than 0, 0 (skip) or 1, l(plot) is 
detected. The tape is stopped for all 
other commands to allow the machine suf- 
ficient time to process the data. When 
using this position, the accuracy of the 
plotted graph is dependent upon the dis- 
tance between points. The MAX POINT 
DISPLACEMENT switch usually will be 
rotated to the FREE RUN position when the 
FREE RUN setting of the TAPE MODE switch 
is used. 

The setting of this switch determines Line 

whether the system functions to cause 
the pen to point-plot or line-plot. The 
POINT position of this switch automatically 
selects the free run filters for high speed 
pen and arm movement. 


This control consists of two concentri- 
cally mounted rotary selector switches. 
The outer switch is set to the length 
which will accommodate the maximum dis- 
tance between two successive points on 
the graph. The inner switch is used in 
conjunction with the outer switch to 
set the time or speed of the plotted 
points. The inner switch can be placed 
to a lower setting than necessary if 
speed rather than maximum accuracy is 
required. 


Outer Switch 
5.0 Inch 
Inner Switch 
Maximum Clock 
wise position 
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CONTROL 


FUNCTION 


REQUIRED 
SWITCH POSITION 


PARITY 

ACTION 


SCALING 

UNITS 


This switch determines the action taken 
by the magnetic tape DATAPLOTTER circuitry 
in the event of a parity error. In the 
STOP position, a parity error switches the 
system to a hold status, the BAD PARITY 
indicator lights, and the fault alarm 
sounds. If a parity error occurs with 
this switch in the REJECT position, the 
accumulators are automatically reset and 
the tape is started forward. Therefore, 
in REJECT, data with parity errors is 
ignored. 

The unit of measure used in plotting is 
set by this switch. The selected posi- 
tion, 1/2 INCH or CM (for centimeter) 
also determines the distance obtained 
for each unit of board offset. If the 
plot is setup in units of half-inches, 
the same values can be used for plotting 
in centimeters by simply changing the 
setting of this switch. No additional 
changes are necessary since all data 
conversion is done in the magnetic tape 
DATAPLOTTER. 


Reject 


1/2 inch 
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APPENDIX A 


1.1 THE EXTENDED FAST FOURIER TRANSFORM 
THE DISCRETE FOURIER TRANSFORM 

Both the Fast Fourier Transform (FFT) and the Extended Fast Fourier 
Transform (EFFT) algorithms compute the discrete Fourier transform or its 
inverse for a given signal. The discrete Fourier transform pair is defined 
by 



N-l 

E 

k=0 


x(k) W' jk 


(A- 1 ) 


N-l 

x ( k ) = E 

j=0 


x(j)Wjj k 


(A-2) 


for j=0, 1, ..., N-l and k=0, 1 N-l. is defined by 

W N = exp(2iri/N) (A-3) 

where i = f-T . 

In the above x(k) represents the time-domain function and X(j) represents 
the frequency-domain function. 


Consider the expression 
/v N - ] ik 

x(j) = E A(k)wf 
k=0 N 


(A-4) 


By comparison with Equations (A-l) and (A-2) 
if A(k) = then x(j) = X*(j) 
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and if 


A(k) = X(j) then x'(j) = x ( k ) 

where the * superscript denotes the complex conjugate. Thus, only Equation 
(A-4) needs to be considered in computing Equations (A-1) and (A-2). 

ALGEBRA OF THE EFFT 

Suppose N, in Equation (A-4), may be written as 
N = LM 

where both M and L are powers of 2. 

The indices j and k in Equation (A-4) may be expressed as 

k = 1 + mL 
j = m' + I'M 

where 1 = 0, 1, - 1;m = 0, 1 M - 1; m* = 0, 1, - 1; 

and T = 0, 1, ..., L - 1. 

In terms of the new indices Equation (A-4) becomes 

x(m' + I'M) = E 1 E 1 A(1 + mL) w (m*+l l M) (1+mL) ( A _ 5 ) 

1=0 m =0 n 

Now 

w (m'+l 'M)(l+mL) = exp pli (m'+TM) (1+mL)] 

- o»n fo-,/m'l mhn , , m \l 

exp n + L + M + ^ m /J 

= ym ' 1 J'l u m ' m 
"N W L W M W 1 
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the term w| m = 1 . 

So Equation (A- 5) becomes 


x(m'+l'M) = y y A(l+mL) W™' 1 W. 1 ' 1 li’ m 

1=0 m=0 IN L M 


(A- 6 ) 


Rewriting Equation (A-6) as 


x(m'+TM) = E 1 Wjj' 1 Wj 1 ' 1 A ( 1 +mL ) wjj' m 

1=0 N L m=0 M 


(A- 7 ) 


In Equation (A- 7) define A-Jm',1) by 

M-l 

E 

m=0 


M-l 

A-| (m ' , 1 ) = E A ( 1 + mL) wJJ 


m 


(A-8) 


Comparing Equations (A-8) and (A-4), then A^m',1) is the discrete 
Fourier transform evaluated at index or frequency m' of the series produced 
by taking every L th sample from A. These samples are produced by beginning 
from the l^* 1 sample of A. There are L of the sequences A^(m',l), each has 
M terms. 

Defining 

A(m',l) by 

A(m',l) = wjj ' 1 A^mM) 

and substituting into Equation (A-7) yields 
^ L- 1 in 

x(m'+1'M) = £ A(m\l) W/ 1 (A-10) 

1 = 0 L 


J 

(A-9) 
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Comparing Equations (A-10) and (A-4) indicates that "x is produced by taking 
the discrete Fourier transform of the M sequences defined by Equation (A-9). 

All FFT transform algorithms are based upon factoring N into products 
and then going through algebra steps as above. The algorithm used as a 
subroutine in the EFFT routine is based upon factoring N into products of 2. 

MATRIX INTERPRETATION 

The algebra of the EFFT has a simple interpretation in terms of matrices 
Consider the original signal term A{1 + mL) in Equation (A-5). The 1 ,m 
indexing is such that A may be stored as a complex matrix with M rows and L 
columns. Thus, m is the row index and 1 is the column index. The discrete 
Fourier transforms defined by Equation (A-8) represent the separate trans- 
forms of each column of the matrix. There are L such columns and each is 
of length M. After transforming, the complex array represented by A^(m',l) 
is stored as a matrix with m' representing the row index and 1 the column 
i ndex . 

The operation defined by Equation (A-9) represents multiplying each 
term of the resulting transforms by an appropriate cos + i sin terms. This 
complex sinusoidal term is called the "twiddle factor." 

The discrete Fourier transform defined by Equation (A-10) is produced 
by transforming separately each row of the complex matrix A. There are M 
rows and each row has L terms. The row index is m 1 and the column index is 1 

DETAILS OF THE EFFT ALGORITHM 

In order to transform extremely long data sequences, the EFFT routine 
that has been implemented uses tapes for storage of the input and output 
signals. Two tapes are required. 

The matrix interpretation of the previous section indicates that the 
general flow of the algorithm must be 

1. The data is input, properly stored as a complex matrix on tape 1. 

2. Transform each column of the matrix on tape 1 and multiply by the 
twiddle factor and store on tape 2. 
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3. Transpose complex matrix on tape 2 and store on tape 1. 

4. Transform each column of the matrix on 1 and store final result 
on tape 2. 

The transpose is required as the algebra indicates than in the first set of 
transforms the columns are operated on while in the second set of transforms 
the rows are operated on. 

The format of the data on the tapes at input or output is as follows: 
Let AR(I,J) and AI(I,J) denote the real and imaginary components of the 
complex matrix respectively. The index I corresponds to the row index where 
1=1, .... NR. The index J corresponds to the column index where 0=1, .... 
NC. Then the matrices AR and AI are stored on the tape as 

Record 1 AR(I,1) I = 1, ..., NR 

Record 2 AI (1 ,1 ) I = 1, .... NR 

Record 3 AR(I,2) 1 = 1 NR 

Record 4 AI (1 ,2) 1 = 1 NR 


Record 2*NC-1 AR(I,NC) I = 1, .... NR 

Record 2*NC AI (I ,NC) 1 = 1, .... NR 

If at input the -row and column dimensions are NR = NR1 and NC = NCI, then 
at output the row and column dimensions are NR = NCI and NC = NR1 . This 
interchange of the dimensions of the AR and AI arrays is due to the trans- 
pose operation. 

The correspondence between the (I,J) indices and the time or frequency 
variables at either input or output is as follows: 

A. If AR and AI represent the real and imaginary parts of a time 
function, then the value of the time variable t is 
t = [( j-1 ) + NC*(I-l)]*At 

for J = 1, NC and I = 1, NR and where At is the time 
separation between points. 
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B. If AR and AI represent the real and imaginary parts of a frequency 
function, then the value of the frequency variable f is 

1. For I = 1, NR/2 

f = [(J-l) + ( I - 1 )*NC]*Af 

for J = 1, NC and where Af is the frequency separation 
between points, i.e., 

a-f - 1 

AT “ NR*NC*At 

2. For I = (NR/2) + 1 , NR 

f = -[(NR*NC- (J-l )- ( I— 1 )*NC]*Af 
for J = 1 , . . . , NC. 
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APPENDIX B 

FILTER TRANSFER FUNCTION IMPLEMENTATION 

Subroutine TRFN computes the filter transfer functions for the 
Butterworth, Chebyshev, and the user specified filters. A detailed 
discussion and derivation of the Butterworth and Chebyshev transfer 
functions are available in Reference 1. 

Butterworth Filter . The transfer function H(s) for the Butterworth • 
filter is implemented in the form 


o n ' 

l+a 1 s+a 0 s <:: + . . . +a s 
1 2 n 

where the a^ 's are constant coefficients built into the software program 
and are given in Table A-l , n is the filter order, and s = jWp (defined 
below). 

The Butterworth filter is defined by the maximally flat response 
function 


|H(J“p) 




f 


i ♦ ^ 


(B-2) 


which has poles given by 


s k = exp[J( ?lc *” - h tf)] . k*l, 2 2n. (B-3) 


The coefficients a^ are derived -by writing Equation (B-l) as 


H(s) = 


(S-S-, ) (s-s 2 ) ( s — s 3 ) (s-s n ) 


(B-4) 
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Table B-l Butterworth Filter Coefficients 


n 

a l 

a 2 

a 3 

a 4 

a 5 

a 6 

1 

1.0 






2 

1 .41421 3562 

2.0 





3 

2.0 

2.0 

1.0 




4 

2.61 3125929 

3.414213562 

2.613125929 

1.0 



5 

3.236067977 

5.236067977 

5.236067977 

3.236067977 

1 .0 


6 

3.863703305 

7.464101615 

9.141620172 

7.464101615 

3.863703305 

1.0 



substituting the poles s ] , s 2> S 3 , . . . s n , performing the indicated 
multiplications in (B-4) to obtain the form of (B-1), and normalizing the 
results such that H(o) = 1. 

Equation B-l is in a frequency normalized form such that 
|H(jo)p)| = 0.707 when = 1 . Thus, for any order Butterworth 
filter, |H(jo)| = 1 and |H{jl)j = 0.707. The normalized form 
may be apdied as a low pass filter by normalizing the frequency 
w to the cut-off frequency w ; i.e., u = 03/03 . Thus, jH(jw)| = 1.0 
when oj = 0 and |H(joj)| = 0.707 when &> = w . The normalized form 
may be applied as a bandpass filter by defining the bandpass center 
frequency u 0 and normalizing u such that ajp = (u-u o )/(cj c -oj o ) for 
o > m and 03 = (ta -o})/ (u 3 -oi ) for u £ w . The upper and lower cut-off 

v U U L U " U 

frequencies are assumed to be symmetrically positioned about a ) Q . 

Summarizing, Butterworth filtering may be obtained by user specifica- 
tion of (1) the filter order, (2) upper cut-off frequency, and (3) center 
frequency of bandpass filter. The center frequency is assumed to be zero 
for a low pass filter. 

Chebyshev Filter . The transfer function H(s) for the Chebyshev 
filter is implemented in the form 


H(s) = 


1 + b,s + boS 2 + . . . + b s n 
i c n 


(B-5) 


where the b. 's are coefficients computed by the software program, 
n is the filter order, and s = jojp (defined below). 

The Chebyshev filter is defined by the equal ripple response 
function 


i H (j“ p )i - i 

^ + C Z n V 


(B-6) 
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where e is the ripple width and C n is the nth order Chebyshev 
polynomial. The pole locations (s k = c k + ju> k ) of Equation (B-6) are 
given by 

\ = ± tanh a s1n[(^l) (f)J (B-7) 

"k = C0S {( 2 H = - L )(f)] (e-8) 

for k = 1, 2, 3, . .., 2n. The parameter a is defined by 

a = 1 sinh" 1 I . ‘ (B-9) 


Implementation of the Chebyshev filter is somewhat more involved than 
than the Butterworth filter because the pole locations (for a given filter 
order) are a function of the input specification e. Thus, note that the 
expression for o k contains the parameter a. Consequently, the coefficients 
b.j must be computed for each problem-oriented application of the filter. 

The equations for this computation are presented in Table B-2. Only odd 
orders are given since only odd order Chebyshev filters are physically 
realizable for equal source and load impedance. 

Equation (B-5) is in a frequency normalized form such that |H(j W p)|ss 
0.707 when w = 1. The normalized frequency, u , is defined 

• r 


ojp = tu/(w c cosh a). 

The normalized form may be applied as a low pass or bandpass filter using 
the frequency normalizing relations (divided by the factor cosh a) de- 
scribed in the discussion of the Butterworth filter. 

Summarizing Chebyshev filtering may be obtained by user specification 
of the (1) filter order, (2) upper cut-off frequency, (3) center frequency 
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of bandpass filter, and (4) ripple width e. 


User Filter . The user may simulate an arbitrary filter by specifying 
selected complex data values of the transfer function and the corresponding 
frequencies. The software program computes the transfer function at de- 
sired frequencies by linearly interpolating the specified data values. 

For example, if the approximation of a transfer function by m straight- 
line segments is desired, the user inputs are m break-point frequencies 
and the associated m complex values of the transfer function. 



k=tanh a=tanh[- sinh' 1 h 
n C J 


Filter Order n Equation 


a, = 1/k 


a, = 


a„ = 


a-> = 


a 2 4b 2 +2a 

(a 2 +b 2 ) 

2a 4 1 
(a 2 4b 2 ) 

(a 2 4b 2 ) 


where a = k/2,b = /3/2 


-> - i j 2a , 2c 

i 1 ~2~j 

a*4b* c^4d 2 


, - 2a , 2c 2.22? 
a ? ~ ~ + ~ 9 — o + a 4b +c 4d 44ac 

a+b c+d (?'+b 2 ) (c 2 +d 2 7“ 


a = - 1 + 1 . 2a42c44ac 

3 a 2 4b 2 c 2 4d 2 (a 2 4b 2 )(c 2 4d 2 


Table B-2. Chebyshev Filter Coefficients 
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Filter Order n Equation 


l+2a+2c 

4 ‘ (aV)(c 2 + d 2 ) 


* 5 ' (.VllcVl 


where a = k cos 36°, b = sin 36 { 
c = k cos 72° , d = sin 72 e 


Table B-2. Chebyshev Filter Coefficients (Continued) 
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